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SU (3) gauge theory coupled to Nf = 2 fermions in the sextet representation is a promising candi- 
date for a technicolor inspired Standard Model extension. In this note the progress in the past few 
years aimed at understanding the non-perturbative properties of the model is reviewed. The main 
difficulties lying ahead in order to make robust conclusions from lattice simulations are outlined. 
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1. Introduction 

Lattice studies of technicolor inspired strongly interacting extensions of the Standard Model 
have seen considerable increase in recent years. In particular SU (3) gauge theory coupled to Nf = 2 
flavors of sextet representation fermions was studied by various methods and techniques. There are 
three main reasons for the interest in this particular model: perturbation theory predicts that it 
lies close to the conformal window and may perhaps give rise to a walking coupling constant, the 
number of would-be Goldstone bosons, if chiral symmetry is spontaneously broken, is exactly three 
and finally the number of flavors is low enough to hopefully give rise to a low S-parameter. Taken 
together these expectations make SU (3) gauge theory with Nf = 2 sextet fermions an appealing 
candidate for strongly interacting extensions of the Standard Model. The role higher dimensional 
representations and the sextet in particular may play in strong dynamics has been emphasized in 
[1,2]. 

The 3-loop /3-function of the model for Nf flavors of Dirac fermions is 



dpi ^ 16?r 2 ^(16tt 2 ) 2 ^(16tt 2 ) 3 

10 „ 250 „ 2857 30475 11425 , 

P l = -11 + -N f , ft = - 102 + —N f , rj 3 = -— + -^Nf - —N} , 

where the first two coefficients [3, 4, 5, 6, 7] are universal and the third [8, 9] is given in MS. 
Results on the 4-loop /3-function with higher dimensional representations can be found in [10]. 
Asymptotic freedom forces Nf < 3.3. The Nf = 1 theory does not have spontaneous symmetry 
breaking because the anomaly breaks the whole flavor group. This leaves the available interesting 
choices as Nf = 2 and Nf = 3 leaving aside half integer Nf values corresponding to Weyl fermions. 

The Nf = 3 model is close to the upper end of the conformal window and is expected to 
be conformal. The fixed point coupling in the 2-loop approximation is gl/An pa 0.085 and it is 
expected that at such a small coupling the 2-loop result is trustworthy. Indeed, taking into account 
the 3-loop correction the fixed point moves to gl/4n pa 0.079 which is only a 7% change. 

The 2-loop result also predicts that the Nf = 2 model is conformal but the fixed point coupling 
gl/4n pa 0.842 is strong enough to doubt the validity of the 2-loop approximation. Again taking 
into account the 3-loop /3-function the fixed point moves to gl/4n pa 0.5 which is a 40% change 
relative to the 2-loop result; it can hardly be considered a small perturbation. In fact the ladder 
resummation in the Schwinger-Dyson approach predicts that the Nf = 2 model is actually chirally 
broken [1, 11]. 

If chiral symmetry is broken for Nf = 2 the expected pattern of chiral symmetry breaking is 
SU (2) x SU (2) — > SU (2) because the sextet representation is complex. This would result in three 
Goldstone bosons precisely the right number for the three massive electro-weak gauge bosons. 

In order for a model to be viable from a phenomenological point of view the ^-parameter 
needs to be not too high otherwise it would conflict with precision data [12]. The phenomenon 
of walking requires the model to be just below the conformal window and achieving this typically 
involves the increase of the fermionic content [13, 14, 15]. For the fundamental representation the 
interesting range is then around Nf = 10 — 14 whereas for the sextet, as we have seen is around 
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Nf = 2. A lower Nf hopefully means a lower S-parameter even with increased dimension for the 
representation, although this statement is also only valid in perturbation theory. 

The above listed appealing features are all based on perturbation theory and/or analogy with 
QCD. Non-perturbative simulations are then needed to clarify which of these expectations do hold. 
Before more detailed questions can be asked the basic property whether in the infrared the model 
is chirally broken or conformal needs to be settled. Ideally, various approaches should yield the 
same conclusion. The various approaches include different physics aspects (thermodynamics, mass 
spectrum, running coupling, etc) and different regularization details (Wilson fermions, staggered 
fermions, improved vs. unimproved, etc). Agreement between different regularizations is of course 
only expected once the results are extrapolated to the continuum. This is especially important 
because in a walking or conformal model the system responds to a change in scale very mildly 
even over many orders of magnitude. If one wishes to resolve this small continuum change of the 
system in a lattice simulation at finite lattice spacing one needs to make sure that discretization 
errors are much smaller than the small expected continuum change. 

In this review three approaches to studying the infrared properties of the Nf = 2 sextet model 
are outlined. At present none of the conclusions from the three approaches are extrapolated reliably 
to the continuum. As we will see currently there are disagreements between the simulation results 
but it might turn out that they are entirely due to discretization errors. 

2. Methods 

There are many observables that are sensitive to the low energy dynamics of the model. The 
physics of chiral symmetry breaking is very different from unbroken conformal symmetry and there 
are many ways to characterize this distinction. 

2.1 Running coupling and step scaling 

The continuum renormalized coupling in the infrared distinguishes directly between the two 
cases of interest: if chiral symmetry breaking takes place the coupling keeps growing towards the 
infrared whereas if conformal and chiral symmetry is intact the coupling approaches a fixed point. 
Different coupling schemes might differ in the value gl of the fixed point but the fact that a fixed 
point exists does not depend on the scheme. Hence any non-perturbatively well-defined scheme 
can be used and any physical scale can play the role of the running scale pL. 

In lattice calculations a popular choice is the Schroedinger functional scheme [16]. In this case 
jU = l/L where L is the finite size of the box. Instead of computing the /3 -function as the response of 
the coupling to an infinitesimal scale change as in the continuum one typically calculates a discrete 
version on the lattice. This discrete /3 -function or step scaling function measures the response of the 
coupling to a finite change, say a factor of 2, in the scale. An infrared fixed point is both a zero of 
the ordinary /3 -function and the discrete step scaling function. As in any lattice calculation a careful 
continuum extrapolation is needed, the step scaling function can be zero at finite lattice spacing but 
this artifact fixed point can disappear for smaller lattice spacing. In this context one lattice spacing 
refers to one lattice size change, say N A — > (2N) 4 , where N = L/a. Changing the bare coupling j8 
is then used to change the renormalized coupling and the continuum limit corresponds to ,/V — > °°. 
The continuum limit is especially important because often times the calculations are pushed to 
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Figure 1: The step scaling function calculated in [17] (left) with thin links indicating an infrared fixed point. 
Using fat links for the fermion action (right) the fixed point disappears [18]. See the text for more details. 

large bare coupling in order to see traces of a possible fixed point exactly where lattice artifacts are 
large. 

The calculation of the running coupling in the Schroedinger functional scheme using Wilson 
fermions was started in [17] for the Nf = 2 sextet model. Using an unimproved (thin link) Wilson 
fermion action a zero of the step scaling function was measured at one lattice spacing corresponding 
to 4 4 — > 8 4 , see left panel of figure 1. Two more lattice spacings corresponding to 6 4 — > 12 4 and 
8 4 — > 16 4 were then added [18] using an improved (fat link) Wilson fermion action, see right panel 
of figure 1. The fixed point disappeared with a possible interpretation that the rougher lattice 
spacing result was an artifact. The gauge action was the same in the two calculations. However 
changing not only the fermion action but the gauge action as well to use fat links resulted in a step 
scaling function with a zero for the lattice spacing corresponding to 6 4 — > 12 4 , see figure 2. A 
possible interpretation is that the absence of the zero previously was the artifact after all [19]. 

Changing the action and/or the lattice spacing led to results so far which show that discretiza- 
tion effects are still there. Clearly a careful continuum extrapolation is necessary with a given 
action in order to decide which finite lattice spacing result is the one prevailing all the way to the 
continuum. A good check of the procedure would be the reproduction of the 2-loop /3 -function for 
small renormalized coupling, carefully extrapolated to the continuum. 

As a cross-check it would be helpful if the running coupling would be calculated in a different 
non-perturbatively well-defined scheme. Reproducing the 2-loop j3 -function for small coupling is 
always a good test for any scheme. For larger coupling two schemes can disagree on the value 
of the coupling but if a fixed point exists for one scheme a fixed point should exist for the other 
scheme too. 

2.2 Thermodynamics 

Another way of addressing the infrared behavior of the model is studying it at finite tempera- 
ture. If chiral symmetry is broken at T = one expects a chiral symmetry restoration temperature 
T c . If the model is conformal in the infrared then as far as chiral symmetry is concerned there 
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Figure 2: The step scaling function from [19] using fat links for the fermion action only (blue) and fat links 
for both the fermion and gauge actions (black). The fixed point is visible again; see the text for more details. 

is no phase transition at all for T > 0. Lattice investigations of thermodynarnical properties are 
complicated by the fact that the lattice system at finite lattice spacing typically has a rich phase 
structure with various types of phase transitions and phase boundaries most of which however hap- 
pens to be regularization specific and as such an artifact with no consequence to the continuum. 
Bulk phase transitions are an example. A careful continuum extrapolation of the findings is hence 
again essential. 

The thermodynamic study of the Nf = 2 sextet model was initiated in [20]. Using unimproved 
rooted staggered fermions in the fixed— N t approach the Polyakov loop, the chiral condensate and 
chiral susceptibility were measured at various quark masses. In the fixed— N t approach one lattice 
spacing corresponds to a given = 1/ (aT), j8 is used to change the temperature and the continuum 
limit is achieved via Nt — >■ °°. A thermal phase transition corresponds to a critical f$ c (N t ) coupling 
for each N t which for large N t scales according to the continuum /3-function; in particular jS c — > °°. 
A bulk phase transition on the other hand is characterized by critical /3 C (A^) couplings which do 
not scale and for large N t approach a fixed value. 

As always with any thermodynamics study finite volume effects need to be under control and 
the quark mass needs to be small enough. Since staggered fermions are used the lattice spacing 
also needs to be small enough in order to avoid dangerous taste violation effects especially because 
the low energy dynamics is very sensitive to the number of massless flavors. 

The critical coupling jS f was determined in [20] from the peak of the chiral susceptibility on 
N t = 4 lattices for two values of the quark mass. The location of the peaks appear to be mass 
independent and is around j8 c « 6.3, see top left panel of figure 3. The N t = 6 result at the same 
two quark masses also from [20] is shown on the top right panel of figure 3. The critical coupling 
moved to j8 c « 6.6. On even finer lattices [21], at N t = 8, the critical coupling moved further, to 
around /3 C = 6.7 with additional small quark masses added, see bottom left panel of figure 3. Again 
the quark mass dependence is quite small. Finally the N t = 12 lattices are preliminary [22] at the 
moment but seem to indicate further increase in j3 c , see the bottom right panel of figure 3. If indeed 
j3 c scales with N t correctly the located phase transitions would correspond to a continuum phase 
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Figure 3: The chiral susceptibility on N t = 4 and N t =6 lattices (top) from [20], and on N t — 8 and N t = 12 
lattices from [21] and [22] respectively (bottom). 
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transition indicating chirally broken symmetry at zero temperature. 

A priori it is not clear how large N t needs to be in order to be in the scaling regime. Most 
importantly the thin link action suffers from possible large taste violation. Unfortunately, these 
effects are not quantified yet. One could in principle reduce them by using smeared actions. In any 
case a continuum extrapolation is necessary. 

2.3 Meson spectrum 

Yet another way of deciding whether the model in the infrared is conformal or chirally broken 
is via studying the mass spectrum as a function of the quark mass. If chiral symmetry is broken 
chiral perturbation theory [23] provides a quantitative description of the quark mass dependence of 
low energy observables. These typically contain polynomial and logarithmic terms. The leading 
formulae are followed by sub-leading terms and it is not clear how small the quark mass needs to be 
in order for the sub-leading terms to be suppressed sufficiently. The order in the chiral expansion 
which needs to be taken into account may also depend on the observable as well. The functional 
form of the leading and sub-leading terms is however known. 

The application of infinite volume chiral perturbation theory requires of course that the finite 
volume of the system is large enough so that increasing it further does not change the measured 
observables. Finite volume effects are exponential and controlled by m % L as usual but how large 
this needs to be is model dependent. Models expected to be close to the conformal window typically 
display much larger finite volume effects than QCD and m K L > 8 — 10 is sometimes necessary. The 
applicability of infinite volume chiral perturbation theory also requires a large f n L value but again 
it is not clear how large it needs to be and the bound probably depends on the observable. Small 
enough lattice spacing is again essential if staggered fermions are used because of taste violation. 

If the model is conformal in the infrared the leading quark mass dependence is power-like and 
fixed by the mass anomalous dimension, M(m) ~ m l '^ l+ ^'. This scaling holds for all masses and 
decay constants as well [24]. For the chiral condensate one has ~ m + cm^^^^ 1 ^. Again it 
is not known a priori how small the quark masses need to be in order to neglect sub-leading terms. 
The functional form of the sub-leading terms is however not known. The scaling formula holds 
in infinite volume and just as with the chirally broken case finite volume effects are needed to be 
small in order to apply them. Finite volume effects are again exponential in m K L where now m K 
simply means the mass of a particle with the same quantum numbers as the pion of QCD. The finite 
size of the system can also be used to study finite size scaling and investigate whether a consistent 
scaling variable can be built out of the quark mass m and L. 

The meson spectrum was first studied in [25] using stout improved rooted staggered fermions 
and Symanzik gauge action. The bare quark mass was varied in the range 0.003 — 0.014, the gauge 
coupling was set to j3 =3.2 while the lattice volume was 24 3 x 48 and 32 3 x 64. The finite volume 
effects were not fully explored in [25] and the measured observables on 32 3 x 64 were interpreted 
as infinite volume results. This assumption was subsequently justified in [26]; see below. 

The 1-loop infinite volume chiral perturbation theory fits for the pion mass, decay constant 
and chiral condensate are shown on figure 4 from [25]. The chiral condensate itself is shown 
and also an independently measured quantity which subtracts the term linear in m using the chiral 
susceptibility. The chiral limit value is of course consistent between the two determinations as it 
should be. The advantage of using the subtracted condensate is that this quantity changes much 
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Figure 4: 1-loop chiral perturbation theory fits to the pion mass (top left), decay constant (top right) and 
chiral condensate (bottom) from [25]. For the pion mass the leading order result is also shown in pink. 



less over the fitted mass range than the condensate itself making the extrapolation to the chiral limit 
more reliable. 

It is clear from figure 4 that both the condensate and the decay constant extrapolate to non-zero 
values in the chiral limit, consistently with spontaneous chiral symmetry breaking in the massless 
theory. 

In order to test the compatibility or incompatibility with a conformal massless theory, leading 
order conformal fits using power-like forms were also performed. The resulting global # 2 /dof 
values for the two fits are 1.24 for the chirally broken hypothesis using 1-loop chiral perturba- 
tion theory and 6.96 for the conformal hypothesis. The global fit in both cases used the same 3 
channels, pion mass, decay constant and chiral condensate. Clearly the analysis in [25] favors the 
chirally broken hypothesis over the conformal one. However it is doubtful whether 1-loop chiral 
perturbation theory is trustworthy in the full fermion mass range for the decay constant; see below. 

As mentioned already the applicability of infinite volume formulae for both the chirally broken 
and conformal scenarios requires large enough lattices. This aspect of the analysis in [25] was 
investigated in [26] where additional simulations at /3 = 3.2 and m = 0.003 were performed on 
48 3 x 96 lattice volumes. At this mass the pion mass, decay constant and the chiral condensate 
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Figure 5: Volume dependence of the pion mass (top left), decay constant (top right) and chiral condensate 
(bottom) for m = 0.003, ]8 = 3.20 and volumes 24 3 x 48, 32 3 x 64 and 48 3 x 96 from [26]. 



are shown on figure 5 together with (an essentially exponential) fit to the volume dependence. As 
can be seen the difference between the 32 3 and 48 3 spatial volumes is very small hence will be 
even smaller for m > 0.003. Consequently in the quark mass dependence fits the 32 3 volumes can 
indeed be used as infinite volume results for all masses. Figure 5 leads to the requirement m n L > 5 
in order to have essentially no further volume dependence (below 1% level). 

Another shortcoming of the analysis in [25] was also addressed in [26]. It is clear from figure 
4 that the effect of the subleading term in the chiral expansion is quite small for the pion mass 
and subtracted chiral condensate. These two quantities are most probably extrapolated to the chiral 
limit reliably. However the decay constant changes quite a bit over the fitted mass range all the way 
to the chiral limit. The difference between the leading order (constant) and subleading (logarithm) 
terms in the chiral expansion is substantial raising the question whether higher order corrections 
need to be taken into account. The data in [25] does not allow fitting 2-loop chiral perturbation 
theory unfortunately. Another option would be to use the 1-loop formulae for a smaller mass range 
where the deviation between the leading and subleading terms is smaller however this would leave 
too few data points for the analysis and would not determine the unknown coefficients with any 
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precision. 

As mentioned already the mass range where the chiral logarithms become significant is a 
priori not known. The polynomial terms in chiral perturbation theory can also be used to fit the 
data, which was performed in [26] together with increasing the statistics at each point. The mass 
range was set to the smaller interval m = 0.003 — 0.008 in order to avoid higher masses where 
the applicability of any small mass expansion is questionable. Since the results from the 48 3 x 96 
volumes at m = 0.003 gave convincing evidence that the 32 3 x 64 results are quite close to the 
infinite volume limit, again the latter volumes were used in the fits for m > 0.003 and the infinite 
volume extrapolated result for m = 0.003; see figure 5. 

In figure 6 polynomial mass dependence is fitted for the pion mass, decay constant, p mass, 
and for what can be called the Higgs mass which is the ++ particle or a-particle (or resonance) 
of QCD. In composite Higgs models, such as perhaps the Nf = 2 sextet model, this particle cor- 
responds to the composite Higgs. This last plot is incomplete in the sense that neither the dis- 
connected propagator nor mixing with glueball states has been measured and is shown only as 
illustration. As shown on figure 6 the fits have acceptable # 2 /dof values. Taken together with the 
non- vanishing values of the decay constant and chiral condensate they indicate that the data is com- 
patible with spontaneous chiral symmetry breaking and hence non-conformal infrared behavior. 

The conformal fit using the same y in the three channels gave even poorer quality fits than 
previously in [25] even with the smaller fitted quark mass range. One might attempt to fit each 
channel separately using separate y anomalous dimensions in the exponent. This results in incom- 
patible values for y, see for instance figure 7 for the tension between the resulting y's from the pion 
mass and decay constant channels. 

Summarizing the results from the meson spectrum one observes that regardless if one uses 
1-loop chiral perturbation theory or polynomial fits a non- vanishing decay constant and chiral con- 
densate is found in the chiral limit with acceptable fit qualities. The power-like conformal fits on 
the other hand are of poor quality or give non-compatible anomalous dimensions if they are not 
constrained to be the same. Hence the analysis in [25] and [26] favor spontaneous chiral symmetry 
breaking over conformal infra red behavior. 

There are several ways to improve on the simulations. At weaker coupling /3 = 3.25 simula- 
tions are ongoing and taste splitting is much improved over j8 =3.2. The conformal analysis could 
be improved by considering subleading terms in the quark mass dependence of the observables. As 
mentioned already finite size scaling could be used combining information from several volumes 
and looking for universal scaling. 

It is worth noting that the small deviation from the leading order behavior m^{m) ~ m hints 
that even if the model is conformal the anomalous dimension is large. From the scaling m K (m) ~ 
m V( 1 +r) one g e ts y~ 1. The result from [19] on the smallness of / is in contradiction with the 
above observation and it remains to be seen how this disagreement gets resolved once the necessary 
improvements are performed to both analyses. 

Though the mass spectrum analysis clearly suggest a chirally broken phase, it is important to 
reduce taste splitting and to perform a continuum extrapolation by using finer lattices. 
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Figure 6: Polynomial fits for the quark mass dependence of the pion mass (top left), decay constant (top 
right), p-mass (bottom left) and the composite Higgs ++ mass (bottom right) from [26]. 



3. Summary and outlook 

The past several years have seen renewed interest in strong dynamics in the context of the 
Standard Model. Part of the reason for it is the realization that higher dimensional representations 
for the techni fermions may solve some of the problems associated with constructions based on the 
fundamental representation. One of the simplest and most appealing constructions is SU (3) gauge 
theory with Nf = 2 fermions in the sextet representation. 

The most basic properties required in order to have any chance for this model to be relevant 
for phenomenology are all non-perturbative. Lattice simulations are an ideal tool to address these 
issues and several methods have been used by several groups. So far there is no final result which 
has been carefully extrapolated to the continuum and passed all necessary consistency checks. The 
three approaches based on the Schroedinger functional, thermodynamics and mass spectrum dis- 
agree in that the thermodynamics and mass spectrum study seems to indicate the model is chirally 
broken whereas the latest incarnation of the Schroedinger functional approach seems to suggest it 
is conformal. 
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Figure 7: Conformal fit to the decay constant (left) and pion mass (right) using separate y anomalous 
dimensions in the exponents; from [26]. As shown the obtained values are not compatible. 

The mass anomalous dimension also presents a disagreement: the Schroedinger functional 
analysis [19] favors a small value y < 0.45 whereas the mass spectrum leads to a much larger 
1 < 7 < 2. Since universality is only expected in the continuum it is entirely possible that the 
disagreements so far are a result of finite cut-off effects and a more systematic approach of the 
continuum is necessary. 

In this review only the Nf = 2 model was mentioned but there are interesting recent quenched 
results [27, 28] and also a thermodynamics analysis of the Nf = 3 case [29]. Both of these cases 
are useful to study because unlike in the Nf = 2 model it is known what the result should be. The 
Nf = model is chirally broken whereas the Nf = 3 model is almost certainly conformal in the 
infrared. Any method used should reproduce these two features. 
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